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Abstract 

This paper describes a measurement of the W boson transverse momentum distribution using 
ATLAS pp collision data from the 2010 run of the LHC at yfs = 7 TeV, corresponding to an 
integrated luminosity of about 31 pb _1 . Events from both W — > eu and W — > \xv are used, 
and the transverse momentum of the W candidates is measured through the energy deposition 
in the calorimeter from the recoil of the W. The resulting distributions are unfolded to obtain 
the normalized differential cross sections as a function of the W boson transverse momentum. 
We present results for p^ < 300 GeV in the electron and muon channels as well as for their 
combination, and compare the combined results to the predictions of perturbative QCD and a 
selection of event generators. 
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I. INTRODUCTION 



At hadron colliders, W and Z bosons are produced with non-zero momentum transverse to 
the beam direction due to parton radiation from the initial state. Measuring the transverse 
momentum (p?) distributions of W and Z bosons at the LHC provides a useful test of 
QCD calculations, because different types of calculations are expected to produce the most 
accurate predictions for the low-pr and high-p^ parts of the spectrum. This measurement 
complements studies which constrain the proton parton distribution functions (PDFs), such 
as the W lepton charge asymmetry in pp collisions because the dynamics which generate 
transverse momentum in the W do not depend strongly on the distribution of the proton 
momentum among the partons. The W pt is reconstructed in W — > Iv events (where I = e or 
\i in this paper). Because of the neutrino in the final state, the W pr must be reconstructed 
through the hadronic recoil, which is the energy observed in the calorimeter excluding the 
lepton signature. This measurement is therefore also complementary to measurements of 
the Z pt, which is measured using Z — > ii events in which the Z pt is reconstructed via the 
momentum of the lepton pair . Although the underlying dynamics being tested are similar, 
the uncertainties on the W and Z measurements are different and mostly uncorrelated. 
The transverse energy resolution of the hadronic recoil is not as good as the resolution on 
the lepton momenta, but approximately 10 times as many candidate events are available 
{{a w ■ BR{W -> iv))l{a z ■ BR{Z ->■ ££)) = 10.840 ± 0.054 [3]). Testing the modeling of 
the hadronic recoil through the W pr distribution is also an important input to precision 
measurements using the W — > li> sample, including especially the W mass measurement. 

In this paper, we describe a measurement of the transverse momentum distribution of W 
bosons using ATLAS data from pp collisions at ^fs = 7 TeV at the LHC [4], corresponding 
to about 31 pb _1 of integrated luminosity. The measurement is performed in both the 
electron and muon channels, and the reconstructed W pt distribution, following background 
subtraction, is unfolded to the true pr distribution. Throughout this paper, p^ is used to 
refer to the reconstructed W Pt and p^ is used to refer to the true W p?- The true W 
Pt may be defined in three ways. The default in this paper is the px that appears in the 
W boson propagator at the Born level, since this definition of p^ is independent of the 
lepton flavor and the electron and muon measurements can be combined. It is also possible 
to define p^ in terms of the true lepton kinematics, with ("dressed") or without ("bare") 
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the inclusion of QED final state radiation (FSR). These define a physical final state more 
readily identified with the detected particles, so we give results for these definitions of p^ 
for the electron and muon channels. For all three definitions of , photons radiated by the 
W via the WW^j triple gauge coupling vertex are treated identically to those radiated by a 
charged lepton. 

The unfolding proceeds in two steps. First, a Bayesian technique is used to unfold the 
reconstructed distribution (p?) to the true distribution {p^) for selected events, taking 
into account bin-to-bin migration effects via a response matrix describing the probabilistic 
mapping from p^ to p^. This step corrects for the hadronic recoil resolution. Second, the 
resulting distribution is divided in each bin by the detection efficiency, defined as the ratio 
of the number of events reconstructed to the number produced in the phase space consistent 
with the event selection. This converts the p? distribution for selected events into the p? 
distribution for all W events produced in the fiducial volume, which is defined by p l T > 
20 GeV, Jt^I < 2.4, p v T > 25 GeV, and transverse mass m T = ^j2p t T p l ^{l — cos(^ — ip v )) > 
40 GeV 5], where the thresholds are defined in terms of the true lepton kinematics. 

The unfolding results in the differential fiducial cross section da^/dp^ , in which the sub- 
script in (Tfid indicates that the cross section measured is the one for events produced within 
the phase space defined above. The electron and muon differential cross sections are com- 
bined into a single measurement via \ 2 minimization, using a covariance matrix describing 
all uncertainties and taking into account the correlations between the measurement channels 
as well as across the p^ bins. The resulting differential cross section is normalized to the 
total measured fiducial cross section, which results in the cancellation of some uncertainties, 
and compared to predictions from different event generators and perturbative QCD (pQCD) 
calculations. 

This paper is organized as follows. Section [XX] reviews the existing calculations and mea- 
surements of p^ . The relevant components of the ATLAS detector are described in Sec- 
tion II I II and the generation of the simulated data used is described in Section IIVI The event 
selection is given in Section |V] and the estimation of the backgrounds remaining after that 
selection is explained in Section |VTJ The unfolding procedure is described in Section IVHI 
Section IVIIII summarizes the systematic uncertainties. The electron and muon channel re- 
sults, the combination procedure, and the combined results are all given in Section HXl We 
conclude with a discussion of the main observations in Section [X] 
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II. QCD PREDICTIONS AND PREVIOUS MEASUREMENTS 



At leading order, the W boson is produced with zero momentum transverse to the beam 
line. Non-zero pr is generated through the emission of partons in the initial state. At 
low pt, this is dominated by multiple soft or almost collinear partons, but at higher px, the 
emission of one or more hard partons becomes the dominant effect. Because of this, different 
calculations of da/dp^ may be better suited for different ranges of p^ . 

At large p^ {p^ ^ 30 GeV), the spectrum is determined primarily by hard parton 
emission, and perturbative QCD (pQCD) calculations at a fixed order of a s are expected 
to predict da/dp^ reliably [6|. The inclusive cross section prediction is finite, but the 
differential cross section diverges as p^ approaches zero. Differential cross sections calculatec 
to 0(a 2 s ) are available for Z/Y production through the FEWZ [J, 8| and DYNNLO 0, 10 1 
programs, and are becoming available for the W. The MCFM generator 11[ can predict 
Px at 0(ot 2 ) through the next-to-leading order (NLO) calculation of the W + 1 parton 
differential cross section. 

As p?f becomes small, contributions at higher powers of a s describing the production of 
soft gluons grow in importance. These terms also contain factors of ln(M^ / '{p^) 2 ) which 



diverge for vanishing p^ . The Px distribution is better modeled in this regime by calcula 



tions that resum logarithmically divergent terms to all orders in a s 



generator 
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13|. The Resbos 



15| resums the leading contributions up to the next-to-next-to-leading loga- 



rithms (NNLL), and matches the resummed calculation to an 0(a s ) calculation, corrected 
to 0(a 2 ) using a /c-factor depending on p T and rapidity, to extend the prediction to large pYf . 
It also includes a non-perturbative parametrization, tuned to Drell-Yan data from several 



experiments 
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161 ]. to model the lowest Px values. 



Parton shower algorithms such as Pythia |17( and Herwig [18( can also provide finite 
predictions of dojdptf in the low-p^ region by describing the soft gluon radiation effects 
through the iterative splitting and radiation of partons. Pythia implements leading-order 
matrix element calculations with a parton shower algorithm that has been tuned to match 
the pf data from the Tevatron [19l-l2ll]. Similarly, the Mc@nlo [22[ and Powheg |23|-|26| 



event generators combine NLO (0(a s )) matrix element calculations with a parton shower 
algorithm to produce differential cross section predictions that are finite for all p^ ■ 
Generators such as Alpgen 



27j | and Sherpa [28] calculate matrix elements for higher 



orders in a s (up to five), but only include the tree-level terms which describe the production 
of hard partons. Parton shower algorithms can be run on the resulting events, with double- 
counting of parton emissions in the phase space overlap between the matrix element and 
parton shower algorithms removed through a veto [^tJ or by reweighting 29, 30 1. Although 
these calculations do not include virtual corrections to the LO process, they are relevant for 
comparison to the highest pt part of the p^ spectrum, which includes contributions from a 
W recoiling against multiple high-pr jets. 

The W pt distribution has been measured most recently at the Tevatron with Run I 
data (pp collisions at y/s = 1.8 TeV) by both CDF [3l| and DO [32| . Both of these results 
are limited by the number of candidate events used (less than 1000), and by the partial 
unfolding which does not take into account bin-to-bin correlations. The present analysis 
uses more than 100,000 candidates per channel and a full unfolding of the hadronic recoil 
which takes into account correlations between bins, resulting in greater precision overall and 
inclusion of higher-p^ events compared to the Tevatron results. 

Although this is the first measurement of the W pt distribution at the LHC, the W — > £u 
sample at y/s = 7 TeV has been studied recently by both the ATLAS and CMS collabo- 
rations. The ATLAS collaboration has measured the inclusive W — > Iv cross section {^J 
and the lepton charge asymmetry in W — >■ events [1.]. The CMS collaboration has also 
measured the inclusive cross section j^], and has measured the polarization of Ws produced 
with pY > 50 GeV, demonstrating that the majority of W bosons produced at large pt in 
pp collisions are left-handed, as predicted by the standard model [34| . 



III. THE ATLAS DETECTOR AND THE pp DATASET 



A. The ATLAS Detector 



The ATLAS detector [35] at the LHC consists of concentric cylindrical layers of inner 
tracking, calorimetry, and outer (muon) tracking, with both the inner and outer tracking 
volumes contained, or partially contained, in the fields of superconducting magnets to enable 
measurement of charged particle momenta. 

The inner detector (ID) allows precision tracking of charged particles within \r)\ ~ 2.5. 
It surrounds the interaction point, inside a superconducting solenoid which produces a 2 T 



axial field. The innermost layers constitute the pixel detector, arranged in three layers, both 
barrel and endcap. The semi-conductor tracker (SCT) is located at intermediate radii in 
the barrel and intermediate z for the endcaps, and consists of four double-sided silicon strip 
layers with the strips offset by a small angle to allow reconstruction of three-dimensional 
space points. The outer layers, the transition radiation tracker (TRT), are straw tubes 
which provide up to 36 additional R — ip position measurements, interleaved with thin layers 
of material which stimulate the production of transition radiation. This radiation is then 
detected as a higher ionization signal in the straw tubes, and exploited to distinguish electron 
from pions. 

The calorimeter separates the inner detector from the muon spectrometer and measures 
particle energies over the range \r)\ < 4.9. The liquid argon (LAr) electromagnetic (EM) 
calorimeter uses a lead absorber in folded layers designed to minimize gaps in coverage. It 
is segmented in depth to enable better particle shower reconstruction. The innermost layer 
( "compartment" ) is instrumented with strips that precisely measure the shower location in 
T). The middle compartment is deep enough to contain most of the electromagnetic shower 
produced by a typical electron or photon. The outermost compartment has the coarsest spa- 
tial resolution and is used to quantify how much of the particle shower has leaked back into 
the hadronic calorimeter. The hadronic calorimeter surrounds the electromagnetic calorime- 
ter and extends the instrumented depth of the calorimeter to fully contain hadronic particle 
showers. Its central part, covering \rj\ < 1.7, is the tile calorimeter, which is constructed of 
alternating layers of steel and scintillating plastic tiles. Starting at \r)\ ~ 1.5 and extending 
to \r)\ ~ 3.2, the hadronic calorimeter is part of the liquid argon calorimeter system, but 
with a geometry different from the EM calorimeter and with copper and tungsten as the 
absorbing material. The forward calorimeters (FCAL), also using liquid argon, extend the 
coverage up to \r)\ ~ 4.9. 

The muon chambers and the superconducting air-core toroid magnets, located beyond 
the calorimeters, constitute the muon spectrometer (MS). Precision tracking in the bending 
plane (R — rf) for both the barrel and the endcaps is performed by means of monitored drift 
tubes (MDTs). Cathode strip chambers (CSCs) provide precision rj — ip space points in the 
innermost layer of the endcap, for 2.0 < \r)\ < 2.7. The muon triggers are implemented via 
resistive plate chambers (RPCs) and thin-gap chambers (TGCs) in the barrel and endcap, 
respectively. In addition to fast reconstruction of three-dimensional space points for muon 
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triggering, these detectors provide y? hit information complementary to the precision r] hits 
from the MDTs for muon reconstruction. 

B. Online Selection 

The online selection of events is based on rapid reconstruction and identification of 
charged leptons, and the requirement of at least one charged lepton candidate observed 
in the event. The trigger system implementing the online selection has three levels: Level 1, 
which is implemented in hardware; Level 2, which runs specialized reconstruction software 
on full-granularity detector information within a spatially limited "Region of Interest" ; and 
the Event Filter, which reconstructs events using algorithms and object definitions nearly 
identical to those used offline. 

In the electron channel, the Level 1 hardware selects events with at least one localized 
region ("cluster") of significant energy deposition in the electromagnetic calorimeter with 
Et > 10 GeV. Level 2 and the Event Filter check for electron candidates in events passing 
the Level 1 selection, and accept events with at least one electron candidate with Et > 
15 GeV. The electron identification includes matching of an inner detector track to the 
electromagnetic cluster and requirements on the cluster shape. The trigger efficiency relative 
to offline electrons as defined below is close to 100% within the statistical uncertainties in 
both data and simulation. 

The online selection of muon events starts from the identification of hit patterns consistent 
with a track in the muon spectrometer at Level 1. For the first half of the data used in this 
analysis, there is no explicit threshold for the transverse momentum at Level 1, but in 
the second half, to cope with increased rates from the higher instantaneous luminosity, a 
threshold of 10 GeV is used. Level 2 and the Event Filter attempt to reconstruct muons in 
events passing the Level 1 trigger using an ID track matched to a track segment in the MS. 
Both apply a pt threshold of 13 GeV for all of the data used in this analysis. The trigger 
efficiency relative to the offline combined muon defined below is a function of the muon pt 
and rj, and varies between 67% and 96%. Due to its larger geometrical coverage, the endcap 
trigger is more efficient than the barrel trigger. The trigger path starting from a Level 1 
trigger with no explicit Pt threshold is slightly more efficient (1-2%) than the one with a 10 
GeV threshold. 
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C. Data Quality Requirements and Integrated Luminosity 



Events used in this analysis were collected during stable beams operation of the LHC in 
2010 at ^/s = 7 TeV with all needed detector components functioning nominally, including 
the inner detector, calorimeter, muon spectrometer, and magnets. The integrated lumin osity 



is 31.4 ± 1.1 pb 1 in the electron channel and 30.2 ± 1.0 pb 1 in the muon channel 36|, 137 ]. 



IV. EVENT SIMULATION 



Simulated data are used to calculate the efficiency for the W — > Iv signal, to estimate the 
number of background events and their distribution in , to construct the response matrix, 
and to compare the resulting normalized differential cross section (1 / ' o^ido^l d-Pr) to a 
variety of predictions. 

The simulated W — > iv events used to calculate the reconstruction efficiency correc- 
tion and to construct the data-driven response matrix are generated using Pythia version 
6.421 (l7j with the MRST 2007 LO* PDF set [38] . The electroweak backgrounds (W -> tv 



and Z/'-f* — > £ + £ ) are estimated using other Pythia samples generated in the same way. 



Simulated tt and sing 



e-top events are generated using Mc@nlo version 3.41 22] and the 
CTEQ6.6 PDF set 39]. For those samples, the Herwig generator version 6.510 18j is used 
for parton showering and Jimmy version 4.1 40j is used to model the underlying event. 
The muon channel multijet background estimate uses a set of Pythia dijet samples with 
a generator-level filter requiring at least one muon with \rj\ < 3.0 and > 8 GeV. The 
multijet background estimate in the electron channel uses a Pythia dijet sample with a 
generator-level filter requiring particles with energy totaling at least 17 GeV in a cone of 
radius AR = \J (Ar/) 2 + (Aip) 2 = 0.05. In both channels, the normalization of the multijet 
background is set by the data. The multijet samples are used to provide an initial esti- 
mate of the background in the electron channel, and to extrapolate data-driven background 
estimates from control data to the signal region in the muon channel. 

In all of the simulated data, QED radiation of photons from charged leptons was modeled 
using PHOTOS version 2.15.4 |4lJ and taus were decayed by TAUOLA version 1.0.2 42]. 
The underlying event and multiple interactions were simulated according to the ATLAS 



MC09 tunes 



43], which take information from the Tevatron into account. Additional in- 



elastic collisions so generated are overlaid on top of the hard-scattering event to simulate 
the effect of multiple interactions per bunch crossing ("pile-up"). The number of additional 
interactions is randomly generated following a Poisson distribution with a mean of two. 
Simulated events are then reweighted so that the distribution of the number of inelastic 
collisions per bunch crossing matches that in the data, which has an average of 1.2 addi- 
tional collisions. The interaction of the generated particles with the ATLAS detector was 



simulated by GEANT4 
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45|. The simulated data are reconstructed and analyzed with 
the same software as the pp collision data. 

The electroweak and top quark background predictions are normalized using the calcu- 
lated production cross sections for those processes. For W and Z backgrounds, the cross 



sections are calculated to next-to-next-to-leading-order (NNLO) using FEWZ |7j, |8|] with the 
MSTW 2008 46] PDFs (see Ref. jjj for details). The ti cross section is calculated at NLO 
with the leading NNLO terms included {^J, setting m t = 172.5 GeV and using the CTEQ6.6 
PDF set. The single-top cross section is calculated using Mc@nlo with m t = 172.5 GeV 
and using the CTEQ6M PDF set. 

We correct simulated events for differences with respect to the data in the lepton recon- 
struction and identification efficiencies as well as in energy (momentum) scale and resolution. 
The efficiencies are determined from selected W and Z events, using the "tag-and-probe" 

n 

method [3J. The resolution and scale corrections are obtained from a fit to the observed Z 
boson line shape. 

Additional W — > tv samples from event generators other than Pythia are used for 
comparison with the measured differential cross section (1/ ' o^ido^j ' dpYf\ The Mc@nlo 
sample used is generated with the same parameters as the ti sample described above. The 
Powheg events are generated using the same CTEQ6.6 PDF set as the main Pythia 
W — >■ lv samples, and Powheg is interfaced to Pythia for parton showering and hadroniza- 



tion. Alpgen version 2.13 [231 matrix element calculations are interfaced to the Herwig 
version 6.510 Q parton shower algorithm, and use Jimmy version 4.31 [40] to model the un 



derlying event contributions. These events are generated using the CTEQ6L1 PDF set 48]. 



Sherpa event generation was done using version 1.3.0 [28], which includes a Catani-Seymour 



subtraction based parton shower model 49], matrix element merging with truncated show- 



ers 



291 ] and high-multiplicity matrix elements generated by Comix 50j. The CTEQ6L1 



PDF set is used, and the renormalization and factorization scales are set dynamically for 



each event according to the default Sherpa prescription. 



V. RECONSTRUCTION AND EVENT SELECTION 

The measurement is performed on a sample of candidate W — )■ Iv events, which are 
reconstructed in the final state with one high-py electron or muon and missing transverse 
energy sufficient to indicate the presence of a neutrino. 

The event selection used in this paper closely follows that used in the inclusive W cross 
section measurement presented in Ref. j^|. The selection in the muon channel is identi- 
cal to that used in the W lepton charge asymmetry measurement in Ref. The event 
reconstruction and W candidate selection are summarized here. 



A. Lepton (e, fx, and v) Reconstruction 

Electrons are reconstructed as inner detector tracks pointing to particle showers re- 
constructed as a cluster of cells with significant energy deposition in the electromagnetic 
calorimeter. This analysis uses electrons with clusters fully contained in either the barrel 
or endcap LAr calorimeter. These requirements translate into \rj e \ < 2.47 with the transi- 
tion region 1.37 < \r) e \ < 1.52 excluded. To reject background (essentially originating from 



hadrons), multiple requirements on track quality and the e 



ectromagnetic shower profile are 
Track quality criteria include 



applied, following the "tight" selection outlined in Ref. J3]. 
a minimum number of hits in the pixel detector, SCT, and TRT, as well as requirements 
on the transverse impact parameter and a minimum number of TRT hits compatible with 
the detection of X-rays generated by the transition radiation from electrons. The energy 
deposition pattern in the calorimeter is characterized by its depth as well as its width in the 
three compartments of the LAr calorimeter, and the parameters are compared with the ex- 
pectation for electrons. The position of the reconstructed cluster is required to be consistent 
with the location at which the extrapolated electron track crosses the most finely-segmented 
part of the calorimeter. Since electron showers are expected to be well contained within the 
LAr calorimeter, electron candidates with significant associated energy deposits in the tile 
calorimeter are discarded. Finally, electron candidates compatible with photon conversions 
are rejected. Although there is no explicit isolation requirement in the electron identifica- 
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tion for this analysis, the criteria selecting a narrow shower shape in the calorimeter provide 
rejection against non-isolated electrons from heavy flavor decays. With these definitions, 
the average electron selection efficiency ranges from 67% in the endcap (1.52 < \rj\ < 2.47) 
to 84% in the central region < 1.37) for simulated W events. 

Muons are reconstructed from tracks in the muon spectrometer joined to tracks in the 
inner detector. The track parameters of the combined muon are the statistical combination of 
the parameters of the MS and ID tracks, where the track parameters are weighted using their 
uncertainties for the combination. Combined muon candidates with \r)\ < 2.4, corresponding 
to the coverage of the RPC and TGC detectors used in the trigger, are used in this analysis. 
To reject backgrounds from meson decays-in-flight and other poorly- reconstructed tracks, 
the pt measured using the MS only must be greater than 10 GeV, and the px measured in 
the MS and ID must be kinematically consistent with each other: 

|p^? s (energy loss corrected) — p T D \ < 0.5 p l ? . (1) 

For both of these requirements, the momentum measured in the muon spectrometer is 
corrected for the ionization energy lost by the muon as it passes through the calorimeter. 
There are no explicit requirements on the number of hits associated with the MS track, but 
the ID track is required to have hits in the pixel detector, the SCT, and the TRT, although 
if the track is outside of the TRT acceptance that requirement is omitted. Finally, to reject 
background from muons associated with hadronic activity, particularly those produced by 
the decay of a hadron containing a bottom or charm quark, the muon is required to be 
isolated. The isolation is defined as the scalar sum of the pt of the ID tracks immediately 
surrounding the muon candidate track (AR < 0.4). The isolation threshold scales with the 
muon candidate px and is ^Pr* < 0-2 Pt- The combined muon reconstruction and selection 
efficiency varies from 90% to 87% as the muon p T increases from 20 GeV to above 80 GeV. 

The transverse momentum of the neutrino produced by the W decay can be approxi- 
mately reconstructed via the transverse momentum imbalance measured in the detector, 
also known as the missing transverse energy (E™ ss ). The E™ 1SS calculation begins from 
the negative of the vector sum over the whole detector of the momenta of clusters in the 
calorimeter. The magnitude and position of the energy deposition determines the momen- 
tum of the cluster. The cluster energy is initially measured at the electromagnetic scale, 
under the assumption that the only energy deposition mechanism is electromagnetic showers 
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such as those produced by electrons and photons. The cluster energies are then corrected 
for the different response of the calorimeter to hadrons relative to electrons and photons, 
for losses due to dead material, and for energy which is not captured by the clustering pro- 
cess. The E™ 1SS used in the electron channel is exactly this calorimeter-based calculation. In 
the muon channel, the E™ 1SS is additionally corrected for the fact that muons, as minimum 
ionizing particles, typically only lose a fraction of their momentum in the calorimeter. For 
isolated muons, the E™ 1SS is corrected by adding the muon momentum as measured with the 
combined ID and MS track to the calorimeter sum, with the calorimeter clusters associated 
with the energy deposition of the muon subtracted to avoid double-counting. In this context, 
muons are considered isolated if the AR to the nearest jet with > 7 GeV is greater than 



5l| and the E T is measured at the 



0.3. Jets are reconstructed using the anti-A^ algorithm 
electromagnetic scale. For non-isolated muons, the muon momentum is measured using only 
the muon spectrometer. In this case, the momentum loss in the calorimeter is kept within 
the calorimeter sum. To summarize, the E^ lss is calculated via the formula 

isolated non— isolated 

*C = -J2 E iv - E - E piv ( 2 ) 

i j k 

In the above, the E x are the individual topological cluster momentum components, exclud- 
ing those clusters associated with any isolated muon, the p> x y are the momenta of isolated 
muons as measured with the combined track, and the p x are the momenta of non-isolated 
muon as measured in the muon spectrometer. In practice, for the electron channel, only the 
first term contributes, but for the muon channel all three terms contribute. 



B. Event Selection 



Candidate W events are selected from the set of events passing a single electron or a 
single muon trigger. Offline, events are first subject to cleaning requirements aimed at 
rejecting events with background from cosmic rays or detector noise. These requirements 
reject a small fraction of the data and are highly efficient for the W signal [3]. Events 
must have a reconstructed primary vertex with at least three tracks with px > 150 MeV. 
They are rejected if they contain a jet with features characteristic of a known non-collision 
localized source of apparent energy deposition, such as electronic noise in the calorimeter. 
Such spurious jets can result in events with large E™ 1SS but which do not contain a neutrino 
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or even necessarily originate from a pp collision. In the electron channel, events are rejected 
if the electron candidate is reconstructed in a region of the calorimeter suffering readout 
problems during the 2010 run 52]. This last requirement results in a ~ 5% efficiency loss. 

After the event cleaning, we select events with at least one electron or muon, as defined 
above, with transverse momentum greater than 20 GeV. In events with more than one such 
lepton, the lepton with largest transverse momentum is assumed to originate from the W 
decay. To provide additional rejection of cosmic rays, muon candidates must point at the 
primary vertex, in the sense that the offset in z along the beam direction between the 
primary vertex and the point where the candidate muon track crosses the beam line must 
be less than 10 mm. 

Finally, we require E™ ss > 25 GeV and transverse mass greater than 40 GeV to ensure 
consistency of the candidate sample with the expected kinematics of W decay. 

After all selections, 112 909 W — > ev candidates and 129 218 W — >■ \w candidates remain 
in the data. The smaller number of candidate events in the electron channel is mostly due 
the lower electron reconstruction and identification efficiency. 



C. Hadronic Recoil Calculation 

The reconstruction of the W boson transverse momentum is based on a slight modification 
of the E™ ss calculation described above. Formally, the px of the W boson is reconstructed 
as the vector sum of the pr of the neutrino and the charged lepton, p^ = p\ + p\. But 
the neutrino pt is reconstructed through the E™ 1SS , and the E™ 1SS is determined in part 
from the lepton momentum, explicitly in the case of W — > fiv events, and implicitly in 
W — > ev events through the sum over calorimeter clusters. Therefore when the pr of 
the charged lepton and E™ 1SS are summed, the charged lepton momentum cancels out and 
the W transverse momentum is measured as the summed pr of the calorimeter clusters, 
excluding those associated with the electron or muon. This part, which consists of the 
energy deposition of jets and softer particles not clustered into jets, is referred to as the 
hadronic recoil R. The reconstructed p^ is denoted p^ and is defined as the magnitude of 
R. 

In this measurement, the exclusion of the lepton from p^ is made explicit by removing all 
clusters with a AR < 0.2 relative to the charged lepton. This procedure leaves no significant 
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lepton flavor dependence in the reconstruction of p§ } so that it is possible to construct a 
combined response matrix describing the mapping from p?f to pfp which can be applied to 
both channels. To compensate for the energy from additional low-pr particles removed along 
with the lepton, the underlying event is sampled on an event-by-event basis using a cone of 
the same size, placed at the same r\ as the lepton. The cone azimuth is randomly chosen 
but required to be away from the lepton and original recoil directions, to ensure that the 
compensating energy is not affected significantly by these components of the event. The 
distance in azimuth to the lepton is required to satisfy A0 > 2 x AR, and the distance to 
the recoil should match A<p > n/3. The transverse momentum measured from calorimeter 
clusters in this cone is rotated to the position of the removed lepton and added to the original 
recoil estimate. Because this procedure is repeated for every event, the energy in the clusters 
in the replacement cone contains an amount of energy from the underlying event and from 
multiple proton-proton collisions ( "pile- up" ) which is correct on average for each event and 
accounts for event-by-event fluctuations. 

VI. BACKGROUND ESTIMATION 

Backgrounds to W — > ev and W — >■ nv events come from other types of electroweak events 
(Z — > ii and W — > rv\ ti and single top events, and from multijet events in which a non- 
prompt lepton is either produced through the decay of a hadron containing a heavy quark 
(b or c), the decay-in-flight of a light meson to a muon, or through a coincidence of hadronic 
signatures that mimics the characteristics of a lepton. Figure [I] shows the expected and 
observed p? distribution in the electron and muon channels, with background contributions 
calculated as described below. 

Electroweak backgrounds (W — > ru, Z — >■ ££, Z — > tt) and top quark production {ti and 
single top) are estimated using the acceptance and efficiency calculated from simulated data, 
corrected for the imperfect detector simulation and normalized using the predicted cross 
sections as described in Section [TV] These backgrounds amount to about 6% of the selected 
events in the electron channel, and to about 10% in the muon channel. The background in the 
muon channel is larger because the smaller geometrical acceptance of the muon spectrometer 
compared to the calorimeter leads to a greater contribution of Z — > fifi events compared to 
Z — )• ee events. Uncertainties on the summed electroweak and top background rates are 6% 
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at low p!p in both channels, rising to 14% above p§ ~ 200 GeV in the muon channel, and 25% 
in the electron channel. The leading uncertainties on these backgrounds at low pt^ are from 
the theoretical model, since the cross sections used to normalize them have uncertainties 
ranging from 4% (for W and Z) to 6% (for tt), and from the PDF uncertainty on the 
acceptances, which is 3% [3j. The integrated luminosity calibration contributes an additional 



3. 



36 



371 ] . Important experimental uncertainties include the energy (momentum) scale 
uncertainty, which contributes about 3% (1%) at low p^ in the electron (muon) channel, 
increasing to about 6% (5%) at high p^. At high (pp > 150 GeV), there are also 
significant contributions for both channels from the statistical uncertainty on the acceptance 
and efficiency calculated from simulated events. 

The multijet backgrounds are determined using data-driven methods. In the electron 
channel, the observed E™ 1SS distribution is interpreted in terms of signal and background 
contributions, using a method based on template fitting. A first template is built from the 
signal as well as electroweak and top backgrounds, using simulated events. The multijet 
background template is built from a background-enriched sample, obtained by applying all 
event selection cuts apart from inverting a subset of the electron identification criteria. The 
multijet background fraction is then determined by a fitting procedure that adjusts the 

normalization of the templates to obtain the best match to the observed E^ 1SS distribution. 

fl 

This method has been described in Ref. [3j, and is applied here bin by bin in pfp. The multijet 
background fraction is 4% at low p^, and rises to ~ 9% at high p^. Uncertainties on this 
method are estimated from the stability of the fit result under different event selections 
used to produce the multijet background templates, by propagating the lepton efficiency 
and momentum scale uncertainties to the signal templates, and by varying the range of the 
E™ 1SS distribution used for the fit. These sources amount to a total relative uncertainty of 
25% at low pP, decrease to 5% at p? ~ 35 GeV, and progressively rise again to 100% at 
high pp, where very few events are available to construct the templates. 

In the muon channel, the multijet background is primarily from semi-leptonic heavy quark 
decays, although there is also a small component from kaon or pion decays- in-flight. The 
estimation of this background component relies on the different efficiencies of the isolation 
requirement for multijet and electroweak events, and is based on the method described in 
Ref. [3]. Muons from electroweak boson decays, including those from top quark decays, 
are mostly isolated, and their isolation efficiency is measured from Z — > fifi events. The 
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efficiency of the isolation requirement on multijet events is measured using a background- 
enriched control sample, which consists of events satisfying all of the signal event selection 
except that the muon transverse momentum range is restricted to 15 < pt^ < 20 GeV and 
the E™ 1SS and my requirements are dropped. The measured efficiency is extrapolated to 
the signal region (p^ > 20 GeV, E™ 1SS > 25 GeV, and my > 40 GeV) using simulated 
multijet events. Knowledge of the isolation efficiency for both components, combined with 
the number of events in the W — >■ \iv candidate sample before and after the isolation 
requirement, allows the extraction of the multijet background. As for the electron channel, 
this method is applied for each bin in p^, with the number of total and isolated candidates, 
as well as the signal and background efficiencies, calculated separately for each bin. The 
isolation efficiency for the background is fitted with an exponential distribution to smooth 
out statistical fluctuations arising from the limited number of events passing all of the event 
selection in the simulated multijet data. 

The multijet background fraction in the muon channel is found to be 1.5% at low 
and decreases to become negligible for pj* > 100 GeV. Uncertainties on the estimated mul- 
tijet background include all statistical uncertainties, including those on both the signal and 
background isolation efficiency measurements. The full range of the simulation-based ex- 
trapolation of the isolation efficiency for the multijet background is taken as a systematic 
uncertainty. Subtraction of residual electroweak events in the control samples is also included 
in the systematic uncertainty but is a sub-dominant contribution. The relative uncertainty 
on the background rate varies between 25% and 80%, with the largest uncertainties for 
p$ < 40 GeV. 

VII. UNFOLDING OF THE p§ DISTRIBUTION 

The unfolding of the pt^ distribution to the p^ distribution is performed in two steps. In 
the first step, the background-subtracted p^ distribution is unfolded to the true p^ distri- 
bution, using the response matrix to model the migration of events among bins caused by 
the finite resolution of the detector. The result of this step is the distribution of dN/dpYp of 
all reconstructed W events. In the second step, this distribution is divided by a reconstruc- 
tion efficiency correction relating the number of reconstructed W events to the number of 
generated fiducial W events within each bin. That correction results in the differential cross 

16 




FIG. 1. Observed and predicted pj, distributions in the electron channel (a) and in the muon 
channel (b). 

section dam/dPr- 



A. Unfolding of the Recoil Distribution 

The response matrix describes the relation between pYj( and p!^, the true and reconstructed 
W pTi respectively. It reflects the physics of the process (hadronic activity from soft and 
hard QCD interactions) as well as the response of the calorimeters to low energy particles. 
This is in principle captured by a response matrix drawn from simulated W — Y iv events, 
but the simulation of both aspects carries significant uncertainty. Therefore, the treatment 
of the response matrix includes corrections from Z data to improve the model. 

The Z — > ee and Z — > fifi data are used as a model for the hadronic recoil response in 
W events because the underlying physics is similar but there are two independent ways to 
measure the pt of the Z, through the hadronic recoil or the px of the charged leptons. The 
lepton energy resolution is sufficiently good that the dilepton px can be used to calibrate 
the hadronic recoil, with the dilepton px standing in for the true px and the hadronic recoil 
remaining the "measured" quantity. One could construct a response matrix purely from 
Z — > ii events, but such a matrix would be limited by the relatively small number of 
Z — )• U events in the 2010 data and residual differences between W and Z kinematics and 
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production mechanisms. To incorporate the best features of both the W simulation and Z 
data models, we introduce a parametrization of the hadronic recoil scale and resolution. Fits 
to the real and simulated Z data using this parametrization are used to correct the simulated 
W response, and the resulting corrected parametrization is used to fill the response matrix 
used for the unfolding. 

Following this logic, the response matrix is built in three steps. A first version of the 
response matrix, denoted M M c, is directly filled from simulated W — > tv signal events as 
the two-dimensional distribution of an d Vr ■ The parametrized response matrix M param 
is also based solely on simulated W — >■ iv events but is constructed from a fit to the recoil 
as described below. The final corrected parametrized response matrix Mp°£ am uses the same 
functional form as M param , but with the fit parameters corrected using the response measured 
in Z — > ££ data. Only Mp°^ m is used in the central value of the measurement, but Mmc 
and M param are used in assessing systematic uncertainties, particularly those arising from 
the response matrix parametrization and the unfolding procedure. 

To facilitate the incorporation of corrections from the Z data, we introduce an analyt- 
ical representation of the detector response to p^ , and approximate Mmc vi & a smearing 
procedure. Decomposing R into its components parallel and perpendicular to the W line 
of flight, R\\ and R±, the response is observed to behave as a Gaussian distribution with 
parameters governed by pYj( and YE T , where T>E T is the scalar sum of the transverse energy 
of all calorimeter clusters in the event. By choosing the coordinate system to align with 
the W line of flight, any scale offset ("bias") is in the i2y direction by construction, and 
the Gaussian resolution function is centered at zero in the R± direction. Specifically, the 
approximated response M param is obtained from the Monte Carlo signal sample as follows: 

Y,Et) = # + G[ b(p™), <j||(#, E£ T ) ], (3) 
R ± (p%,EEr) = G[ 0, a ± ( p w,YE T ) ]. (4) 

where G denotes a Gaussian random number, and its parameters b, an and a± are the 
Gaussian mean and resolution parameters determined from fits to the simulation. The bias 
is described according to b(p^) = b + bi \[p^ , independently of EEr- The resolutions follow 
<7||(#,E£ T ) = <7||,o(#) + <7||,i(Pt) x v/E^ and a ± (pf,YE T ) = a ±fi (pf) + a ±A (pf) x 
y/EEr, where the pYp dependence indicates that the fit is performed separately in three 
regions of p^ (pf < 8 GeV, 8 < p^ < 23 GeV, and pf > 23 GeV). The separation of the 
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fit into regions of p^ improves the quality of the fit. 

With the parametrization defined, it is possible to build up a response matrix from a 
set of events using a smearing procedure. Given the p^ and HEt of each event, R\\ and 
R± can be constructed using random numbers distributed according to Equations [3] and EH 
Then pfp ls reconstructed from Rn and R±, and the results are used to fill the relationship 
between p^, and p^. Applying this procedure to the simulated signal sample results in the 
approximate response matrix M param . 

Corrections to this parametrization are derived from Z — > 11 events by applying the same 
procedure to both real and simulated Z events and using the measured decay lepton pair 
momentum p{^ as the estimator of the true Z boson transverse momentum. The hadronic 
recoil calculated as described in Section IV CI has no dependence on the lepton flavor, and 
consistent response is observed in Z — > ee and Z — > events. Therefore we fit the 
combined data from both channels to minimize the statistical uncertainty. The corrected 
smearing parameters are defined as follows: 



^W^corr 


_ ^VK,MC _|_ /£«,data _ 


-f- MC ), 


(5) 


W'.corr 
Ci, 


= cr,, + (cr,. 


— cr,,' J, and 


(6) 


W.corr 


W,MC . / «,data 

= a ± ' +{a ± ' 




(7) 



Above, fr w ' data and b ee,MC are determined as a function of p^, and then used as a function 
of p?f\ ftWicorr anc | ^vk,mc are f unc tions of pYjt throughout. All resolution parameters are 
functions of the reconstruction-level Y>Et- This defines the final, corrected response matrix 
^pamm use d in the hadronic recoil unfolding. 

The parametrization of the bias and resolution parameters in W and Z simulation are 
illustrated in Figures E^a)-|4(a). For these, the bias and resolution are defined with respect 
to the true (propagator) W and Z momenta. The simulated and data-driven bias and 
resolution parameters in Z events are displayed in Figs. [2(b) -IU(b). For these, the bias and 
resolution are defined with respect to the reconstructed dilepton p T . In Figs. [2(a) and[2j(b), 
the bias parametrization is shown only over the range which determines the fit parameters, 
but the parametrization describes the data well up to p^=300 GeV. 

The response matrix is constructed using the following bin edges, expressed in GeV: 

• Reconstruction-level distribution: 0, 4, 8, 15, 23, 30, 38, 46, 55, 65, 75, 85, 95, 107, 
120, 132, 145, 160, 175, 192, 210, 250, 300. 
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FIG. 2. (a) Parametrization of the recoil bias as a function of the vector boson transverse momen- 
tum, b(p^ ,Z ), in W simulation (open squares, solid line) and Z simulation (solid circles, dashed 
line), (b) Parametrization of the recoil bias as a function of the reconstructed lepton pair trans- 
verse momentum, b(p^f), in Z simulation (dashed line) and data (solid squares, shaded band). The 
shaded band shows the uncertainty on the fit. 

• Unfolded distribution: 0, 8, 23, 38, 55, 75, 95, 120, 145, 175, 210, 300. 



The reconstruction-level binning enables more detailed comparisons between data and sim- 
ulation before unfolding, and allows a more precise background subtraction as a function of 
Pt- It has been used in Figure [TJ The bin edges at the unfolded level provide a purity of 
at least 65% across the p^ spectrum, which is large enough to ensure the stability of the 
unfolding procedure. The bins are still small enough to keep the model dependence of the 
result, which enters through the assumption of a particular p?£ shape within each bin, to 
a sub-leading contribution to the overall uncertainty (see the description of the systematic 
uncertainties in Section IVlIIj) . The purity is defined as the fraction of events where the event 
falls in the same bin when the bin edges are defined using p^ as it does when the bin edges 
are defined using p% . 

The unfolding of the hadronic recoil is performed by means of the iterative Bayesian 
algorithm (53)], where the p?f distribution predicted by the simulation is used as first as- 
sumption of the true p^ spectrum, and iteratively updated using the observed distribution. 
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FIG. 3. (a) Parametrization of the recoil resolution a\\(EEr) in W simulation (open squares, solid 
line) and Z simulation (solid circles, dashed line), (b) Parametrization of the recoil resolution 
<7||(££t) in Z simulation (dashed line) and data (solid squares, shaded band). The shaded band 
shows the uncertainty on the fit. 
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FIG. 4. (a) Parametrization of the recoil resolution <7j_(EEt) m W simulation (open squares, solid 
line) and Z simulation (solid circles, dashed line), (b) Parametrization of the recoil resolution 
(j±^{YiEt) in Z simulation (dashed line) and data (solid squares, shaded band). The shaded band 
shows the uncertainty on the fit. 
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This procedure converges after three iterations. 

The statistical uncertainty on the unfolded spectrum is obtained by generating random 
replicas of the reconstruction- level data. First, the distribution from simulation is scaled 
to have an integral equal to the number of events observed in data. For each trial, the number 
of events in each bin is fluctuated according to a Poisson distribution with a mean set by 
the original bin content. The unfolding procedure is used on the fluctuated distribution, 
and the p^ distribution from the same set of simulated events is subtracted from the result. 
The resulting ensemble of offsets is used to fill a covariance matrix describing the impact of 
statistical fluctuations on the result, including correlations between the bins introduced by 
the unfolding procedure. 

Systematic uncertainties receive contributions from the quality of the response parametriza- 
tion approximation, i.e. from the difference between Mmc and M param ; from the statistical 
precision of the data driven corrections defining Mp°£ am ; and from the unfolding procedure 
itself. Their estimation is described in Section IVIII1 

B. Efficiency Correction 

The W — > iv candidate event reconstruction efficiency is subsequently unfolded by divid- 
ing the number of events in each bin of p^ by the detection efficiency correction factor for 
that bin. The correction factor accounts for trigger and detection efficiencies, as well as the 
migration of events in and out of the acceptance due to charged lepton and E!f lss resolution 
effects. It is defined as the ratio of the number of reconstructed events passing all selection 
in each bin to the number of events produced within the fiducial volume in that same bin. 
Note that any migration between bins has already been accounted for by the hadronic recoil 
response unfolding. The efficiency correction is based on the ratio calculated from simulated 
W events, and corrected for observed differences between simulated and real data in the 
trigger and reconstruction efficiencies as well as in the lepton momentum and resolution (see 
Section ITVl) . The corrections for discrepancies between data and simulation are applied as a 
function of the reconstructed lepton kinematics in each bin of pYf . The fiducial volume in the 
denominator is defined by the truth-level kinematic requirements p? > 20 GeV, |?^| < 2.4, 

> 25 GeV, and uit > 40 GeV. For the default, propagator- level p^ measurement, the 
lepton kinematics and transverse mass are defined at the QED Born level, i.e., before any 
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final state QED radiation. For the "dressed" lepton version of the measurement, the charged 
lepton momentum is the sum of its momentum after all QED FSR and the momenta of all 
photons radiated within a cone of AR = 0.2 around the lepton. The cone size is chosen to 
match the cone size used for the lepton removal in the definition of R. The "bare" lepton 
version uses only the charged lepton momentum after all QED FSR. 

In the electron channel, the efficiency rises from ~ 60% at low pf to ~ 80% at pf ~ 
100 GeV, and falls towards ~ 70% at the upper end of the spectrum. In the muon channel, 
the efficiency rises from ~ 80% to ~ 90%, then falls to ~ 80% in the same p^ ranges. 

The efficiency correction carries systematic uncertainties induced by the imperfect mod- 
eling of the lepton trigger and reconstruction efficiencies, by the acceptance of the E™ 1SS cut, 
and by the finite statistics and physics assumptions of the signal simulation sample. Their 
estimation is described in Section IVIIII 

VIII. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties arise from the background subtraction procedure, from the recoil 
response model and unfolding procedure, and from lepton reconstruction and calibration 
uncertainties. Theoretical uncertainties also enter, to a lesser extent. Different strategies are 
used for the various uncertainties according to the nature of the uncertainty and whether 
it is introduced before, during, or after the hadronic recoil unfolding. Accordingly, the 
uncertainties are evaluated by using an ensemble of inputs with the nominal response matrix, 
an ensemble of response matrices with the nominal input, or by simple error propagation, 
respectively. The uncertainties on this measurement are represented as covariance matrices, 
so that correlations between the bins can be included. 

A. Background Subtraction Uncertainties 

The systematic uncertainties associated to the background subtraction are estimated 
by generating an ensemble of pseudo-experiments in which the background estimates have 
been fluctuated within their uncertainties. The full analysis chain is repeated for each 
pseudo-experiment and the spread of the unfolded results defines the associated uncertainty. 
Electroweak, top, and QCD multijet contributions are treated separately, except that the 
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luminosity uncertainty is treated as correlated between the electroweak and top backgrounds. 
Background subtraction is performed before the unfolding, and the unfolding redistributes 
the background among the pYjt bins, so the covariance matrices representing the uncertainties 
on the backgrounds have nonzero off-diagonal elements. 

The electroweak and top backgrounds contribute 0.6% (0.4%) to the measurement un- 
certainty at low pj? in the electron (muon) channel, and up to 4% at high pYjt in both 
channels. The multijet background in the electron channel contributes ~ 0.5% uncertainty 
for pif < 50 GeV, which gradually rises to 4% at p^ ~ 200 GeV, eventually contributing 
15% in the highest p^ bin. In the muon channel, the multijet background induced un- 
certainty has a maximum of 2% at p^ ~ 30 GeV, which corresponds to the peak of the 
background rate, and contributes ~ 0.6% on average in the rest of the spectrum. 

B. Hadronic Recoil Unfolding Uncertainties 

Systematic uncertainties associated to the response matrix are classified in two categories. 
In the first category, the impact of a given source of uncertainty is estimated by comparing 
the unfolded distribution obtained with the nominal response matrix, to the result obtained 
with a response matrix reflecting the variation of this source. The statistical component of 
the difference is assessed by varying a given input of the response matrix construction to 
generate a set of related variations of the response matrices. Repeating the analysis with 
these leads to a set of varied unfolded results, and the induced bias is averaged in each bin 
of the p^ distribution. The associated systematic uncertainty is defined from the spread 
of the distribution of the results, and is taken as a constant percentage across all p^ bins, 
represented as a diagonal covariance matrix. 

By comparing results obtained from the initial Monte Carlo response matrix Mmc with 
results obtained from the parametrized response matrix M param , the response parametriza- 
tion is found to induce an uncertainty of 2.4% in the electron channel and 2.0% in the 
muon channel. The input generator bias is estimated by reweighting the true p^ distribu- 
tion given by the Pythia sample to the Resbos prediction, generating the corresponding 
response matrix and comparing the result to the nominal result, leading to a systematic 
uncertainty of 1.2% in the electron channel, and 0.9% in the muon channel. Note that the 
starting assumption for the Bayesian unfolding is simultaneously modified in the same way, 
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so that this uncertainty includes both the effect of modifying the distribution underlying the 
response matrix and the assumption of a prior for the unfolding. In addition, it was verified 
that reweighting the input pYjt assumption according to the actual measurement result and 
repeating the procedure does not affect the result beyond the uncertainties quoted above. 
Lepton momentum scale uncertainties also enter through the Z-based recoil response correc- 
tions, because is used in place of the true pf , but this amounts to less than 0.2% in both 
channels. As described above, these numbers are taken constant across the p^ spectrum. 

The second category deals with the uncertainties associated to the data-driven corrections 
to the response parametrization. In this case, we generate an ensemble of random correction 
parameters by sampling from the distribution defined by the statistical uncertainties on the 
central value of the parameters returned by the fit. For each parameter set the correspond- 
ing response matrix is generated. The treatment is then the same as for the background 
uncertainties: the analysis chain is repeated for each configuration, and the spread of the 
unfolded bin contents defines the associated uncertainty in each bin. 

In this category, the data driven correction to the recoil bias and resolution induces an 
uncertainty of ~ 1.6% for p^ < 8 GeV, has a local maximum of ~ 2.6% at p^ = 30 GeV, 
and contributes less than 1% in the remaining part of the spectrum. The uncertainty related 
to the TjEt rescaling is 0.2% at low p? , rising to 1% at the high end of the spectrum. These 
numbers are valid for both channels, as the data driven corrections are determined from 
combined Z — > ee and Z — >• samples, as described in Section IVH Al 

Finally, the bias from the unfolding itself is found by folding the pjf distribution of 
simulated W — > tv events passing the reconstruction-level selection using Mmc and then 
unfolding it using the same response matrix. The original p^ distribution is subtracted from 
the unfolded one, and the size of the bias relative to the original distribution is taken as the 
systematic uncertainty from the unfolding procedure. The folded distribution is used for 
Pt instead of the found p^ distribution to avoid double-counting the statistical uncertainty. 
The resulting uncertainty is less than 0.5% in all bins, except for the highest-p^ bin in the 
electron channel, where it is 1%. 
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C. Efficiency Correction Uncertainties 



In the electron channel, the main contributions to the acceptance correction uncertainty 
are the reconstruction and identification efficiency uncertainty, and the electron energy scale 
and resolution uncertainties. The identification efficiency contributes 3% to the measurement 
uncertainty across the p? spectrum. The scale and resolution uncertainties contribute 0.5% 
at low , rising to 10% at pY? ~ 100 GeV, and decreasing to 6% at the high end of the 
spectrum. 

In the muon channel, the trigger efficiency uncertainty contributes 1% across the spec- 
trum. The reconstruction efficiency contributes 0.7% at low , linearly rising to 2% at 
pYX ~ 300 GeV. The scale and resolution uncertainties contribute 0.5% at low p? , rising to 
2% at pf ~ 120 GeV, and decreasing to 1% towards pYf, ~ 300 GeV. 

The uncertainty associated to the recoil component of E™ ss (the first term of Eq. [5J mi- 
nus any clusters associated with an electron) is estimated as above, by generating random 
ensembles of resolution correction parameters within the precision of the Z-based calibra- 
tion. For each parameter set in the ensemble, the E™ ss distribution is re-generated and the 
corresponding efficiency correction is re-calculated. The width of the resulting distribution 
of efficiency corrections is taken as the uncertainty. This source contributes less than 0.3% 
across the pY/ spectrum in both channels. 

In both channels, the Monte Carlo statistical precision is 0.5% at low pY/ and rises to 
4% towards pY/ ~ 300 GeV. The generator dependence of the efficiency is estimated by 
comparing the central values found for Pythia and Mc@nlo, and found to be smaller 
than 0.2%, apart from the last bin where it reaches 1%. Finally, following Ref. |2|, the PDF 
induced uncertainty on the efficiency correction is at the level of 0.1% and neglected in this 
analysis. 

IX. RESULTS 

A. Electron and Muon Channel Results 

The efficiency-corrected distributions resulting from the two unfolding steps are normal- 
ized to unity, and the bin contents are divided by the bin width. In the normalization step, 
uncertainties that are completely correlated across all of the bins, such as the uncertainty 
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FIG. 5. Electron and muon fiducial differential cross sections as a function of . The error bars 
include all statistical and systematic uncertainties except the 3.4% uncertainty on the integrated 
luminosity, which is common to the two measurements and cancels in the ratio. 



on the integrated luminosity, cancel. The resulting normalized differential fiducial cross 
section, {1/ a^Aj^da^/ dp^) is given in Table [I] for both the electron and muon channels, 
together with the statistical and systematic uncertainties. The differential cross section is 
calculated with respect to three definitions of p^f and the fiducial volume, corresponding to 
different definitions of the true lepton kinematics: the first uses the Born level kinematics, 
the second uses the "dressed" lepton kinematics calculated from the sum of the post-FSR 
lepton momentum and the momenta of all photons radiated within a cone of AR = 0.2, and 
the third ( "bare" ) uses the lepton kinematics after all QED radiation. 

Instead of normalizing the efficiency-corrected distributions to unit integral, they can also 
be divided by the integrated luminosity of the corresponding data to yield the differential 
fiducial cross section da^/dp^ . The resulting differential fiducial cross sections, with the 
fiducial volume defined by the Born-level kinematics, are shown in Fig. [5j Error bars include 
both statistical and systematic uncertainties, but not the uncertainty on the integrated 
luminosity, which is common to both measurements. 
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TABLE I. The normalized, differential cross section (1 / 'ofidXcfofid/ 'dpS ') , measured in W — > ev and 
W —?• [iv events, for different definitions of pj, . The Born-level definition ("propag."), the analysis 
baseline, ignores the leptons and takes the W momentum from the propagator. The "dressed" and 
"bare" definitions of pjf are calculated using the momenta of the leptons from the W decay. In the 
"dressed" case, the charged lepton momentum includes the momenta of photons radiated within a 
cone of AR = 0.2 centered around the lepton. In the "bare" case, the charged lepton momentum 
after all QED radiation is used. The factor p is the power of ten to be multiplied by each of the 
three cross section numbers for each channel. It has been factorized out for legibility. 
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B. Combination Procedure 

After correcting the electron and muon distributions to the common fiducial volume 
using the efficiency corrections described in Section IVHl we combine the resulting differen- 
tial fiducial cross sections da^/dp^ by x 2 minimization. The combination is based on the 
distributions with defined by the W propagator momentum because QED final state 
radiation causes differences between the electron and muon momenta that makes a consis- 
tent combination based on other definitions unfeasible. To build the \ 2 i the uncertainties 
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on the two measurements are sorted according to whether they are correlated between the 
two channels or not, and a joint covariance matrix describing the uncertainty on both mea- 
surements is constructed. Using this covariance matrix, we define a x 2 between the two 
measurements and a common underlying distribution. This \ 2 is minimized to find the 
combined measurement, which is the best estimate of the common underlying distribution. 
Specifically, the x 2 to be minimized is defined as: 



where X is the vector of 2N elements containing the two iV-bin distributions to be com- 
bined, concatenated: X = {Xf, X^\ X% , ...,X%}. The vector X = X n ; Xi, X n } 
contains two copies of the combined measurement {Xi}. The joint covariance matrix C is 
described in the next p arag raph. The \ 2 minimization is performed analytically, following 
the prescription in Ref. [541 ] . yielding the {^}. 

The joint covariance matrix C has IN x IX elements and is constructed from four sub- 
matrices: 



The N x N covariance matrices C e and are the covariance matrices for the electron and 
muon measurements, respectively, and contain all sources of uncertainty on the measure- 
ments. The off-diagonal blocks C e/t are identical and reflect the sources of uncertainty that 
are correlated between the channels. 

The 2N x 2N covariance matrix is constructed from the two N x N matrices for each 
source of uncertainty individually, and the resulting set of 2N x 2N matrices is summed. For 
sources of uncertainty uncorrelated between the channels, the 2N x 2N covariance matrix 
is constructed by copying the N x N matrices to the corresponding diagonal blocks C e and 
(jn p or uncertainties that are correlated between the channels, the diagonal blocks are still 
filled by copying the covariance matrices from the individual channels. The off-diagonal 
blocks are filled using the assumption that the channels are 100% correlated, so that the 
correlations between bins are identical for both channels. That determines the correlation 
matrix, which sets the magnitudes of the covariance matrix entries relative to the magnitude 
of the diagonal entries. The diagonal entries, which are the squares on the uncertainties on 
each bin, are taken as the geometrical average of the values for the two channels. 



X 2 = (X-X) T C- 1 (X-X) 



(8) 




(9) 
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The statistical uncertainties on the unfolded measurements are uncorrelated because the 
W — > ev and W — >■ \iv candidate data samples are statistically independent. The system- 
atic uncertainties induced by the subtraction of the estimated background are uncorrelated 
between the channels, except for the uncertainties on the luminosity and predicted cross 
sections used to normalize the electroweak and top quark backgrounds. Because the same 
hadronic recoil response matrix is used for both channels, the uncertainties associated with 
it are fully correlated between the channels, except for the small contribution from the lep- 
ton momentum resolution. The efficiency corrections for each channel are independent, so 
the associated uncertainties are uncorrelated between the channels. 



C. Combined Results and Comparison with Predictions 



The x 2 minimization yields a x 2 /d.o.f. of 13.0/13, demonstrating good agreement between 
the electron and muon results. The combined differential cross section, normalized to unity, 
is shown compared to the prediction from RESBOS in Fig. EJ The RESBOS prediction, which 
combines resummed and fixed order pQCD calculations, is based on the CTEQ6.6 PDF 
set J39 1 and a renormalization and factorization scale of mw- Resbos performs the fixed- 
order calculation at NLO (0(a s )), and corrects the prediction to NNLO 



[0(a 2 )) using a 

fc-factor calculated as a function of the boson mass, rapidity, and pt |13l-ll5f|. Table HT1 gives 
the same information numerically, including the separate contribution of different classes of 
uncertainty. 

In Figured the combined result (1/ 'ofidX^fid/ 'dPr) is compared to a selection of predic- 
tions from both pQCD and event generators. The DYNNLO predictions are from version 1.1 



of the program (9 , Ll0[] - The prediction from the MCFM program is produced as a calculation 
of da^/dp^ for W + 1 parton events and uses MCFM version 5.8 11]. The leading order 
calculation for W + 1 parton production is 0{a 3 ) and the NLO calculation is 0(a 2 s ), so the 
predictions are comparable to other 0(a s ) and 0(a 2 ) predictions of p^ for p^ > 5 GeV, the 
minimum jet px threshold in the calculation. Both of the pQCD calculations are normalized 
by dividing the prediction in each bin by the inclusive cross section prediction calculated in 
the same configuration as the differential cross section, and both have the renormalization 
and factorization scales set to mw- The 0(a s ) predictions use the MSTW2008 NLO PDF 
sets, and the 0(a 2 s ) predictions use the NNLO MSTW2008 PDF set |4fj|. The uncertainty 



30 




FIG. 6. Normalized differential cross section obtained from the combined electron and muon 
measurements, compared to the Resbos prediction. 

on the pQCD predictions comes mostly from the renormalization and factorization scale 
dependence, and studies indicate that it is comparable in magnitude to the 10% and 8% 
observed for p% predictions at 0(a s ) and 0(a 2 ) in Ref. j^j]. 

The DYNNLO and MCFM predictions do not include resummation effects and are not 
expected to predict the data well at low pYp because of the diverging prediction for vanishing 
Pr . Therefore, the lowest bin (p^ < 8 GeV) is omitted from Fig. [7J The two programs 
predict similar distributions at the same order of a s . The 0(a s ) prediction from both 
calculations for the fraction of the distribution above p? ~ 23 GeV is about 30% too low 
on average, similar to the NLO event generators. The 0(a s ) prediction from FEWZ [7, 
is not shown in Fig. [7J but is in agreement with those from DYNNLO and MCFM. The 
discrepancy between the predictions and the measurement appears when normalizing to the 
inclusive cross section and would be compensated by a large but unphysical contribution in 
the first bin. The ratio moves closer to unity in the high pY}( range. The 0(a 2 s ) predictions 
agree better with the data than those at 0(a s ). They are within 15% of the data for all p? . 

The predictions of the event generators Pythia, Powheg, Alpgen, Sherpa, and 
Mc@nlo are based on the simulated samples described in Section IIV1 Since Powheg and 
Alpgen can be interfaced with more than one parton shower implementation, the nota- 
tions Powheg+Pythia and Alpgen+Herwig are used to make the choice explicit. The 
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TABLE II. Measured pip using combined electron and muon data, with all uncertainties shown by 
source. 



pYp Bin (l/crfld)(dcrfl ( j/dpy ) ResponseMatrix Backgrounds Efficiency Statistical Total 
[GeV] (GeV^ 1 ) uncert. (%) uncert. (%) uncert. (%) uncert. (%) uncert. (%) 



— 8 


5.510 


i n— 2 
10 


1.91 


0.26 


0.76 


0.22 


O AO 

2.48 


8 — 23 


2.512 


10~ 2 


1.69 


0.28 


0.87 


0.24 


2.42 


23 — 38 


6.766 


10~ 3 


3.20 


0.57 


1.28 


0.57 


4.31 


38-55 


2.523 


10~ 3 


2.34 


0.65 


1.44 


0.84 


3.78 


55-75 


1.025 


10~ 3 


1.78 


0.74 


1.74 


1.19 


4.09 


75-95 


4.263 


10~ 4 


1.61 


1.15 


2.13 


1.91 


4.94 


95 - 120 


1.896 


10~ 4 


1.98 


1.94 


2.67 


2.68 


5.99 


120 - 145 


7.985 


10~ 5 


2.84 


3.30 


3.16 


4.78 


7.91 


145 - 175 


3.710 


10~ 5 


1.98 


2.66 


3.66 


5.72 


9.31 


175-210 


1.692 


10" 5 


2.00 


3.72 


3.84 


7.75 


10.56 


210-300 


4.803 


10~ 6 


2.69 


7.81 


4.26 


9.28 


14.40 



Pythia, Resbos, Sherpa, and Alpgen+Herwig predictions describe the measurement 
within 20% over the entire range. For pY^ < 38 GeV, the data indicate a softer spectrum 
than these predictions. For 38 < p^ < 120 GeV, the data distribution exceeds the Resbos 
prediction and undershoots the Sherpa prediction, but agrees with the Alpgen+Herwig 
and, to a lesser extent, pure Pythia predictions. For p? > 120 GeV, Pythia and Resbos 
agree in predicting a softer spectrum than Alpgen+Herwig and Sherpa, but the data 
provide no significant discrimination among these predictions. 

Powheg+Pythia and Mc@nlo, the NLO event generators interfaced with parton 
shower algorithms, provide a reasonable description of the data for p^ < 38 GeV, but 
both underestimate the data starting at p^ w 38 GeV, with a deficit gradually increasing 
to nearly 40% at high pY^ . 

Finally, we compare the combined result to the measurement of (l/afid)(<^fid/^pf ) de- 



• 6. 



scribed in Ref. [2|. The W and Z have different masses and couple differently to quarks, 
so the results cannot be directly compared, but the ratios of the measured to predicted dis- 
tributions for a common model can be used to qualitatively assess the agreement between 
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FIG. 7. Ratio of the combined measurement and various predictions to the Resbos prediction for 
(l/cr M )(da M /dpf), using (a) the 0(a s ) and 0(a 2 s ) predictions from DYNNLO and MCFM, and 
using (b) the predictions from Alpgen+Herwig, Mc@nlo, Powheg+Pythia, Pythia, and 
Sherpa. The statistical uncertainties on the generator distributions are negligible compared to 
the uncertainty on the measurement and are not shown. 

the two measurements. The ratios of the W and Z distributions in data to their respective 
Resbos predictions are overlaid in Fig. |HJ In spite of the different techniques and uncer- 
tainties characterizing both measurements, the ratios display similar trends as a function of 
p^, the true boson p T . 



X. CONCLUSIONS 

The W transverse momentum differential cross section has been measured for p^ < 300 GeV 
in W — > iv events reconstructed in the electron and muon channels using the ATLAS de- 
tector. The W — >■ iv candidate events are selected from pp collision data produced at 
a/s = 7 TeV, corresponding to approximately 31 pb -1 from the 2010 run of the LHC. 

The measurement is compared to a selection of predictions. The Alpgen+Herwig, 
Pythia, Resbos, and Sherpa predictions match the data within 20% over the entire 
p^ range. Mc@nlo provides the closest description of the data for p^ < 38 GeV, but 
Mc@nlo and Powheg+Pythia both underestimate the data at higher p% . Fixed-order 
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FIG. 8. The ratio of {1/ a^){dcr^/ dp^) as measured in the combined electron and muon data 
to the Resbos prediction, overlaid with the ratio of {1/ a^)(da^/ dp^ ) measured as described in 
Ref. [2|] to the Resbos prediction. 



pQCD predictions from the DYNNLO and MCFM programs agree very well with each 
other. They predict fewer events at high p^ at 0(a s ) but the agreement with the measured 
distribution is significantly improved by the 0{a 2 s ) calculations. 

A comparison of the W and Z data relative to the prediction from a given theoreti- 
cal framework displays similar features across the measured transverse momentum range, 
supporting the expected universality of strong interaction effects in W and Z production. 

Although the measurement is limited by systematic uncertainties over most of the spec- 
trum, the dominant uncertainty sources can be constrained with more integrated luminosity. 
With the integrated luminosity available from the 2011 run now in progress, future measure- 
ments should be able to measure da^/dp^ to at least double the current range in pYp . With 
improved statistical and systematic uncertainties, it should also be possible to measure the 
ratios of the W to Z and W + to W~ differential cross sections as functions of the boson p?, 
which will further test the predictions of QCD. 
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